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Abstract 



I We show that a baryon asymmetry can be generated by dissipative effects during warm inflation via a supersym- 
T— I -metric two-stage mechanism, where the inflaton is coupled to heavy mediator fields that then decay into light species 
[through B- and CP- violating interactions. In contrast with thermal GUT baryogenesis models, the temperature during 



dnflation is always below the heavy mass threshold, simultaneously suppressing thermal and quantum corrections to the 



^~>]inflaton potential and the production of dangerous GUT relics. This naturally gives a small baryon asymmetry close to 



■the observed value, although parametrically larger values may be diluted after inflation along with any gravitino over- 
|abundance. Furthermore, this process yields baryon isocurvature perturbations within the range of future experiments, 
■making this an attractive and testable model of GUT baryogenesis. 
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The inflationary paradigm [l[ has been extremely suc- 
cessful in explaining the flatness and homogeneity of the 
observable universe, as well as providing an origin for Cos- 
mic Microwave Background (CMB) anisotropics and the 
Large Scale Structure (LSS). While models with a single 
slow-rolling scalar field have been extensively explored in 
the literature, these haven proven hard to embed within 
UV-completions of the Standard Model (SM) such as su- 
pergravity /string theory constructions, which typically suf- 
fer from the so-called 'eta-problem' thatprecludes sufii- 
ciently long periods of inflation (see e.g. [2[). 

Warm inflation [1,13 (see also (Hj) offers an attractive 
solution to this problem by taking into account dissipative 
effects, which not only damp the infiaton's motion and al- 
low for longer periods of slow-roll but also act as a source 
of light particles that may lead to a 'graceful exit' into 
a radiation-dominated era. Furthermore, when the tem- 
perature of the radiation exceeds the Hubble rate during 
inflation, T > H, thermal inflaton fluctuations become 
the main source of density perturbations, typically yield- 
ing lower inflationary scales than cold scenarios. 

From a supersymmetric (SUSY) two-stage mechanism 
0, @] , in which interactions between the inflaton and the 
light particles are mediated by heavy 'catalyzer' flelds X, 
whose masses are above the temperature of the radiation 
bath, T < mx/100, the flatness of the inflaton potential 
is safe from both quantum and thermal corrections. As 
supersymmetry is, however, inefficient in cancelling time 
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non-local processes, this allows for strong dissipative ef- 
fects with moderately large field multiplicities [J] . 

Dissipation is naturally an out-of-equilibrium process, 
as annihilation of the resulting particles cannot efficiently 
're-populate' the classical background condensate. It is 
then natural to ask whether a baryon asymmetry may 
be produced during inffation through dissipative effects, 
by incorporating B- and CP-violating interactions in the 
two-stage mechanism described above, thus satisfying the 
Sakharov conditions Q. The structure of the two-stage 
interactions then suggests a parallel with thermal baryo- 
genesis models in GUT scenarios, in which the B- and 
CP-violating decays of heavy GUT bosons may occur in 
an out-of-equilibrium fashion once the temperature of the 
universe drops below their mass threshold [8|. 

In this Letter, we show that heavy GUT states can me- 
diate dissipative processes during warm inflation leading to 
a baryon asymmetry (see e.g. [9| for related studies). An 
attractive feature of this mechanism is the generation of 
an asymmetry at T <C Mqut-, avoiding the production 
of undesirable relics, like magnetic monopoles, that are 
generic in thermal scenarios. As we show below, the result- 
ing asymmetry is naturally small in the low-temperature 
regime, despite the 0(1) couplings typically required for 
successful warm inflation, as opposed to both thermal and 
non-thermal GUT baryo genesis scenarios where these 
are necessarily suppressed. 

We also show that the resulting baryon-to-entropy ratio 
exhibits thermal fluctuations that lead to baryon isocurva- 
ture perturbations in the CMB spectrum and that ma y b e 
accessible with upcoming experiments such as Planck 11| , 
a feature which is lacking in many of the proposed models 
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of baryogenesis. Warm baryogenesis thus constitutes the 
first example of a consistent warm inflation scenario where 
a baryon asymmetry is produced and a viable model of 
GUT baryogenesis. 

We start by considering a GUT-hke extension of the 
SUSY two-stage mechanism of that includes interac- 
tions between the inflaton, described by the scalar compo- 
nent of a superfield and a set of mediator superfields 
Xa, a = I,--' ,Nx- These are coupled to light degrees 
of freedom (dof) described by a set of Nf superfields Qi 
and Li, where the former carry a non-zero baryon number 
(or B L charge) bi and the latter correspond to non- 
baryonic species. For simplicity, we will loosely refer to 
these as 'quark' and 'lepton' superfields, although our dis- 
cussion applies to more generic scenarios. The relevant 
superpotential can be written as: 



W = 



(1) 



where a sum over the heavy and light field indices is im- 
plicit. As generic in SUSY GUT models, the form of the 
Yukawa terms in Eq. ([T]) is such that there is no consistent 
assignment of baryon number to any of the Xa fields, which 
may then decay in a B- violating fashion. The couplings h^^ 
and A^^ are complex and Nf > 3 in order to ensure viola- 
tion of C and CP, whereas the ga couplings are real. The 
mass of the mediator supermultiplet components during 
inflation is given by = 2g1^tp^ , where (0) = ip/ \/2, with 
a negligible splitting between the fermionic and bosonic 
components due to soft SUSY breaking during inflation. 
We will also assume the masses of the Qi and Li compo- 
nents are well below the temperature during warm infla- 
tion and may thus be discarded. 

As shown in the literature (see e.g. Q), the main pro- 
cesses contributing to dissipative effects in warm inflation 
correspond to interactions between the scalar components 
in a superpotential of the form Eq. (U]), so that the dom- 
inant source of a baryon asymmetry will be the produc- 
tion of squark fields. The relevant terms are included in 
the scalar potential Vs = '^glWWa? + 'i9a<p'^ xlQ^i ^3 + 



\]iq\lj) + h.c, where Xa, and Ij denote the scalar com- 
ponents of the Xa, Qi and Lj superfields, respectively. A 
baryon asymmetry will then arise from the difference be- 
tween the particle production rates for scalar quarks and 
antiquarks. In [l^], these were computed in the adiabatic 
regime from the thermal Wightman self-energy for each 
light particle species i, yielding for the time derivative of 
the energy density. 



.(d) 

PI ~ 



dt'- 



-iujp{t — t') 



-^■2l{p,t,t') 



where p and uip are the momentum and energy of the light 
fields. The difference between this quantity for baryonic 
and anti-baryonic species then sets the rate at which a net 
baryon asymmetry is produced, and can be compard to 



the total particle production rate in order to determine 
the baryon-to-entropy ratio produced by warm inflation. 
Here we will outline the main steps of this computation, 
describing the more technical details in the appendix for 
the interested reader. 

Firstly, one should note that in the low-temperature 
regime, where thermal corrections to the inflaton poten- 
tial are suppressed, the main contribution to dissipative 
particle production comes from one-loop effects involving 
virtual heavy Xa bosons, since production of on-shell states 
is Boltzmann suppressed for rria ^ T. This is inherently 
different from most baryogenesis models in the literature 
involving the decay of on-shell heavy states, where the 
out-of-equilibrium condition results from the inbalance be- 
tween direct decay and inverse decay processes fisj . 

The leading particle production process is in this case 
where ai_j = qi_j,lij, via virtual x field 
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pairs. Since the inflaton field is varying (slowly), the re- 
sulting particles cannot efficiently annihilate and fully give 
away their energy back into the background condensate. 
This leads to a net particle production and makes dissi- 
pative effects inherently out-of-cquilibrium, as signaled by 
the non-zero value of the right-hand side of Eq. (O for 
1^ 7^ (see appendix). 

Although the leading contributions to the light parti- 
cle self-energy arise at one-loop order, the associated di- 
agrams correspond to 'squared' tree-level diagrams, and 
hence cannot contribute to the baryon asymmetry accord- 
ing to the theorem by Nanopoulos and Weinberg The 
diagrams contributing to the squark self-energy up to 2- 
loop order arc then illustrated in Fig. 1, with analogous 
diagrams for the anti-squark self-energy including the com- 
plex conjugate couplings. We have computed these contri- 
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Figure 1: Feynman diagrams contributing to the qi self-energy up to 
2-loop order. The filled circles correspond to couplings Ig^^p^ with 
the inflaton background field. 



butions in the low-temperature regime, T <C m^, justifying 
a low-momentum approximation for the heavy boson prop- 
Note that the Xa propagators are 'dressed'. 
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i.e. involve a resummation over light field loop corrections. 
These correspond to the non-vanishing decay width of the 
heavy bosons, ' 'v^hich induces particle production de- 
spite their virtual nature in the low-temperature regime 
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(see appendix) . We then obtain for the difference between 
the squark and anti-squark self-energies: 



Nx 



ASfi = 8 ^ / d^fc famllablm{TT[XlXlhi,hl] 



(3) 



where k denotes the momentum of the Xa propagators in 
the main loop, fa{p,k,t,t') includes the main loop fac- 
tors and is explicitly given in the appendix, while Iab{p, k) 
is the imaginary part of the 'triangle' loop integrals. In 
the low-temperature regime for the heavy particles, a low- 
momentum approximation applies and we obtain 



1 



m^Iab{p,k) ~ - — ( 1 - - 



Ip^ -{p- kf 



(4) 



From Eqs. ^ and (|3]), it is easy to conclude that, as in 
thermal GUT baryogenesis models 0], one needs at least 
two distinct heavy fields, with either different masses or 
different Yukawa couplings, in order to produce a non- 
vanishing baryon asymmetry. In particular, note that the 
trace Tr[A^A^hbht] is real for a = & and that the lab and 
Iba contributions cancel each other unless rUa 7^ Tit. 

The total particle production rate can be computed by 
adding the self-energy of all the light baryonic and non- 
baryonic fields, for which it is sufficient to consider the 
lowest-order diagrams in Fig. 1. This gives: 



d'k fa 



Tr[h,ht]+Tr[A:Ar]+Tr[A„At] 



(5) 



where the first two terms correspond to the self-energy 
of the light baryonic species and the last term to the non- 
baryonic fields. The internal and external momentum inte- 
grations, as well as the time integration in Eq. ©, can then 
be performed following the procedure in [l2[ . The momen- 
tum dependence of the triangle loop integrals lab yields an 
additional factor T^/ml with respect to the leading order 
result, while changing the overall value by a factor ~ 3.5, 
which we obtained numerically (see e.g. \XM)- This cor- 



responds to an inherent momentum cut-off in Eq. @ due 
to Boltzmann suppression of the heavy mediators in the 
low-T regime, such that p, < T. 

For concreteness, let us consider a case with Nx — 1 
degenerate fields of mass mx and an additional multi- 
plet of mass mx + Amx, with Amx/mx ^ 1. We also 
assume without loss of generality that all Qi superfields 
carry the same baryonic charge, b, and that the Yukawa 
couplings only differ by a phase, i.e. ha = j/ilVe*^" and 
Xa = |A|J7e*"", with U and V denoting unitary matrices. 
We find the asymmetry ratio to be 



■id) 

Pb 

■id) 
PR 



3.5 
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Nx TTij^ mx 



, (6) 



where S denotes an effective CP- violating phase. Note that 
this is independent of the number of light species Nf and 
only mildly dependent on Nx- 



During inflation, the radiation density is determined by 
both the dissipative effects and the Hubble friction, with 
evolution equation pr + 411 = TifP'. The dissipative 
coefficient, computed in the low-T regime [iMIiSjHIj is of 
the form T = C^T'^/ip^, where ~ 0.16|/i|47VxA^J,,,y 
for \h\ ~ |A|, with N decay denoting the number of available 
decay channels for the heavy mediators. Similarly, a frac- 
tion rT sources a net baryon number density ub ~ rib — n^. 
Since pn = {g,,Tr^ /30)T'^, with denoting the number of 
relativistic dof, the entropy and baryon number densities 
are given by: 



s + 3Hs 
u'b + iHuB 



45C(3) QB T(^2 

T 
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(7) 



with g B giving the number of baryonic dof. Note that this 
includes fermionic fields, which are thermally produced 
despite the dominant scalar nature of the dissipative pro- 
14| . In the slow-roll regime, these quantities quickly 



cesses 



approach a steady state solution, yielding for the baryon- 
to-entropy ratio rjs — hb/s, 



rjs w3.5 



45C(3) gBb sin (5 \h\^\X\^ Nx - 1 T"^ Amx 



g, |/i|2 + 2|A|2 Nx m\ mx 



(8) 

In warm inflation models, radiation typically comes to 
dominate the energy density Q, ending inflation during 
the slow-roll regime, such that Eq. ^ gives the observed 
value if no significant entropy production or washout pro- 
cesses occur in the post-inflationary evolution. In partic- 
ular, one can ensure that electroweak sphaleron processes 
(l5| conserve the produced asymmetry by generating a net 
B — L number density during warm inflation. Also, if 
the 'catalyzer' fields are heavy (Majorana) right-handed 
neutrinos with L- violating interactions [l6[, an inflation- 
ary lepton asymmetry may later be converted into baryon 
number by sphaleron processes. 

Apart from numerical factors, the main difference be- 
tween the baryon asymmetry produced by dissipation and 
the corresponding result for out-of-equilibrium decay of 
heavy bosons resides in the /m\ suppression in the low- 
T regime, so that the heavy fields are not excited by back- 
ground dissipation. This is an appealing feature, since it 
not only prevents the generation of GUT relics that could 
overdose the universe but also leads to a small baryon- 
to-entropy ratio, at the same time ensuring that thermal 
corrections are under control. Taking for instance b = 1/3, 
\h\ ~ |A|, Nx > 1 and the MSSM values for 53 and g*, 
we then obtain 



ry, « 8.9 x 10-"|/i| 



T/mx 
0.01 



Am J 



sin (5 
0.015 } V 0.025 



(9) 



Note that the moderately large field multiplicity typically 
required in warm inflation models will not affect the ra- 
tio gBjg* significantly. The best present estimate for the 
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baryonic asymmetry comes from Big-Bang nucleosynthe- 



sis 



nl, 7.2x10-" < 77s < 9.2x10-", at 95% C.L. Our re- 



sult ^ may thus be naturally within the observed window 
despite the 0(1) couplings typically required in warm in- 
flation scenarios, with mass degeneracy and a CP-violating 
phase at only the few percent level. On the other hand, 
thermal GUT baryogenesis scenarios require strongly sup- 
pressed couplings in order to reproduce the observed asym- 
metry. Similarly, in non-thermal scenarios [lo| . the para- 
metric resonance that excites the heavy bosons is efficient 
only for sufficiently long lifetimes. 

In some warm inflationary scenarios, radiation is al- 
ways subdominant even though dissipation sustains the 
required number of e- folds |J]. One then expects addi- 
tional entropy production through conventional reheating 
processes, as well as possibly from the late decay of any 
light moduli present during inflation. Notice, however, 
that moduli dynamics may be modified due to the ther- 
mal radiation bath with T > H. Although this may dilute 
any undesired gravitino abundance |18j . it will also dilute 
the asymmetry produced by dissipation, but this may be 
compensated with larger CP phases, couplings and non- 
degencracies in the heavy mediator spectrum. 

Given that mx oc cfi, Eq. also implies 77^ oc 
so that thermal fiuctuations of the infiaton field will be im- 
printed on the baryon-to-entropy ratio. Although baryons 
are subdominant during inflation, they become a signifl- 
cant component of the energy density at late times and 
will contribute to CMB anisotropics and LSS. Although 
such baryon isocurvature perturbations (BIP) also arise 
in other scenarios [l9| . they are in this case fully (anti- 
)correlated with adiabatic perturbations, since both origi- 
nate from inflaton fluctuations. These can be obtained by 
perturbing the inflaton equation, (jj + {3H + T)(f> + = 
0, and are coupled to the temperature fluctuations via 
the evolution equation for p^. BIP are conventionally 
measured by the ratio Bb ~ Sb/Ci where C = —HSp/p 
is the gauge- invariant curvature perturbation and Sb = 
Spb/pb - iS/4)SpR/pR = 5-qs/ris In the slow-roll 

regime, on supcrhorizon scales, we obtain after some alge- 
bra: 



(i + g)2(i + 7g) ' (''^^ 

where we define Q — T/3H, ~ {m'\,/2){V^/VY , rj^ = 
m'^pV^^/V and = m^pV^ / {V (j)) ■ This thus yields an 
additional observable that may be used to probe the con- 
sistency of warm inflation models, along with the spectral 
index, the tensor-to-scalar ratio and non-gaussianity pa- 
rameters. In particular, we obtain Q 



Bf 



775-1 



3e,^— 2i)^+3a-^ 
3Q 3t^+r-i^-6a^ 



Q < 1 

g > 1 



(11) 



so that Bb is generically at most of the same order of mag- 
nitude as the deviations from scale invariance, with ns = 



0.968 ±0.012 (68% C.L.) [2l|, being further suppressed for 
strong dissipation. For example, for a quartic potential 
one flnds Bb — —0.096, — 0.007/g in the weak and strong 
dissipation regimes, respectively. The most recent WMAP 
analysis of cold dark matter anti-correlated isocurvature 
perturbations [2l| . taking into account that Vlc/^b — 5, 
yields \Bb\ < 0.34 (95% C.L.), for Bb < 0, while accord- 
ing to an earlier analysis —0.53 < Bb < 0.43 (95% C.L.) 
[22|. This is thus generically consistent with the expected 
amount of BIP, even for g ^ 1, while Planck should im- 
prove these bounds by an order of magnitude [llj. 

Generating a baryon asymmetry through dissipative ef- 
fects during warm inflation thus exhibits several attractive 
features. By keeping the temperature during inflation be- 
low the GUT scale, it provides an appealing model of GUT 
baryogenesis, with no unwanted relics and a small baryon- 
to-entropy ratio with no unnaturally suppressed couplings 
and CP-violating angles, at the same time keeping thermal 
corrections to the inflaton potential under control. Fur- 
thermore, a superpotential of the form in Eq. ([T]) is nat- 
ural in SUSY GUT constructions where the inflaton lies 
within a diagonal subgroup of an adjoint Higgs representa- 
tion and the heavy mediators correspond to its off-diagonal 
components, as in the D-brane model considered in (23j . 
These constructions are also natural arenas for warm infla- 
tion due to the large held multiplicity available, ensuring 
sufficiently strong friction effects to sustain the required 
number of e-folds. 

Although in this work we have focused on the baryon 
asymmetry generated by the varying inflaton fleld, it is 
possible in some models of warm inflation for excitations 
above the background condensate, ~ (0) + (fii, to also 
play a signiflcant role. The slow roll conditions require 
ry^ < 1 + g, where Q ~ T/3H, which imply m^/H < 
\/l + g. Since T > H during warm inflation, one could 
have 7770 < T if g ^ 1 in the early stages of inflation, such 
that the process x ~^ ^'i4'i contributes to P^ and hence 
the (pi particle states become a non-negligible component 
of the radiation bath. They are, however, unstable and 
decay via 0i — ?> crfcr'^ through virtual x fields. If dissipa- 
tion becomes strong towards the end of inflation, g ^ 1, 
as is the case for typical inflationary potentials [J], it is 
possible to have 777.0 > T, depending on the evolution of 
Lp and T, so that the (j)i states decay out-of-cquilibrium 
and produce a net baryon number. Given the common 
origin and structure of the relevant interactions, the con- 
tributions from background dissipation and (pi decay to 
the baryon asymmetry will necessarily have the same sign 
and similar magnitude, so that our earlier results are ro- 
bust despite the dependency of the latter contribution on 
the details of the inflationary dynamics. 

One of the crucial aspects of the model presented in this 
Letter is the fact that dissipation provides a novel mech- 
anism for satisfying the Sakharov out-of-equilibrium con- 
dition 0. Although we have examined a nearly thermal 
equilibrium regime, warm inflation is general to any non- 
equilibrium inflationary scenario governed by fluctuation- 
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dissipation dynamics Q. Furthermore, dissipative pro- 
cesses may also play an important role in the post- in- 
flationary universe, e.g. during cosmological phase transi- 
tions, so that we may extrapolate from our analysis a more 
general mechanism of dissipative baryogenesis. 

Warm baryogenesis thus provides an unprecedented link 
between inflationary physics and the observed baryon asym- 
metry, a connection that may be tested in the near future 
and possibly shed a new light on two of the most important 
problems in modern cosmology. 
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Appendix A. Particle production rates 

The scalar interactions involving the background infla- 
ton field. If, the mediators Xa and the light fields ct^ = qi, h 
relevant for the dissipative dynamics are of the form (omit- 
ting the indices for simplicity): 



(A.1) 



where M = \/2gLp. Following [l2|, the particle production 
rate of a-particles in the radiation bath due to interactions 
with the x-fields is given by: 



,{p) = Im 



g-lLj(p)(t-t') 

dt' E2i(p,i,t') 



(A.2) 



with the Wightman self-energy given, to leading order, by: 



(p,t,i') 



l&g^h^M'^ 



d^k 



dujj 



{2n 



-iuj-2 {t—t') 



xG,,2i(p -k,t- t')Gl^{k,ui)^\uji - us) 
xG^_Jk,w3)v''(^3-^2)Gx,6i(k,c^i) , (A.3) 

where G^^abO^yUj), a,b = 1,2, are the fuUly dressed fi- 
nite temperature Schwinger-Keldysh x-propagators. In 
the low-temperature regime, T <|; m^, the main contri- 
butions to this integral arise from virtual low-momentum 
modes, so that we can use the approximate forms: 



G^.a(k,c.) 



2r;,(k,c.) 



1 + n{uj)] 



(A.4) 



where n(uj) is the Bose-Einstein distribution function and 
is the finite temperature width of the x-ficld, which has 
been computed in [IJ] at one-loop order. In the slow-roll 
regime one can approximate the Fourier transform of the 



background inflaton field (p'^ito) ~ 2i(p(p2TTd' {uj) . Also, the 
light CT-propagators can be Fourier transformed to yield: 



Ga,2l{P - k, w) 



-i[l + n{Lu)] 



27T 



2uj° 



(A.5) 



where sgn(a;) is the sign function. Substituting these re- 
sults into Eq. (jA.3[) . we get: 



S,^2i(p,i,i') = J d^'k f{p,k,t,t') 



(A.6) 



where 



f{p,k,t,t') 



GAig'^ f duj2 1 



(27r)3 ml J 2^ 

X {T^{k,L02)[l+n{L02)]}" 

<)(*-*')n(^-_j}(A7) 



+ e 



corresponds to the function appearing in Eqs. and ([S]). 
To obtain the net particle production rate it remains to 
perform the time integration in Eq. (|A.2p . for which we 
may use the following identity: 



7r<5(w2+^^±<_k) + 



.(A.8 



This yields: 



na{p) 



2,;,2 



^32g^h^M^^ 



d^k 



ml J {27t)^uj-u; 



p'^P-k 

{r^(k, -u; u;_^) [I + n{-u; u;_^)] }" 

X [1 + 7l(aj^_k)] + • 

+ {r^(k, -cj; + c.^_k) [1 + n{-u; + u;_^)] }' 



(A.9) 



The total particle production rate in an expanding uni- 
verse Per + 4:Hpc = T(^^ can then be obtained by integrat- 
ing over the momenta of the produced cr-particles: 



T 



1 



d^p 
(2^ 



'^p-riaip) 



(AlO) 



These results can then be used to compute the production 
rate of each species. For the difference between squark 
and antisquark production rates one needs to include the 
additional momentum dependence and the associated cou- 
pling structure arising from the triangle loop diagrams in 
the above expressions, according to Eq. (|4]), while for the 
total particle production rate it suffices to add the one-loop 
contributions of all species, with the relevant couplings. 
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